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ABSTRACT 

We argue that the recent groundbreaking discovery by Badenes et al. (2009) of a nearby {D w 50 pc) white 
dwarf-neutron star (or black hole) binary (SDSS 1257+5428) with a merger timescale < 500 Myr implies that 
such systems are common; we estimate that there are of order ~ 10^ in the Galaxy. Although subject to 
large uncertainties, the nominal derived merger rate is Fmw ^ 5 x lO'^^yr"' in the Milky Way, just ^3-6 
and ~ 20-40 times less than the Type la and core-collapse supernova (SN) rates, respectively. This implies 
that the merger rate is ^ 0.5- 1 x lO'^Gpc"^ yr"' in the local universe, ^ 5000- 10000 times more than the 
observed (beaming-uncorrected) long-duration gamma-ray burst (GRB) rate. We estimate the lower limit on 
the rate in the Galaxy to be Fmw ^ 2.5 x 10~^yr~' at 95% confidence. We briefly discuss the implications 
of this finding for the census of long- and short-duration GRBs and their progenitors, the frequency of tight 
binary companions to Type-Ib/c SN progenitors, the origin of ultra-high energy cosmic rays (UHECRs), the 
formation of rapidly rotating neutron stars and ^2-3 M0 black holes, the census of faint la-like SNe, as well as 
for upcoming and current transient surveys (e.g., LOSS, PTF, LSST), and for high- (LIGO) and low-frequency 
(LISA) gravitational wave searches. 
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1. INTRODUCTION: THE ARGUMENT 

The recent paper by Badenes et al. (2009) describes their 
discovery as part of the SWARMS survey of a probable white 
dwarf-neutron star/black hole (WD-NS/BH)^ binary (SDSS 
1257+5428) at a distance of D « 48^|"pc, with orbital period 
of K, 4.6 hr, and radial velocity semi-amplitude of w 320km 
s"'. The WD is of spectral type DA, has a g-band magni- 
tude of 16.8, and has a cooling age of fcooi ~ 2.0 ± l.OGyr. 
The system masses are Mwd — 0.92;!;o3| and MNs/BHsin(0 ~ 
1.62;|;q25 M0. The transverse velocity of the system is w 
1 1 km s"' and the total spatial velocity is plausibly ~ 120km 
s"' . For / = 60°, the semi-major axis is k, 0.01 AU. The orbital 
period and masses imply a merger timescale for the system of 
+342 jy[yj. ^jj^ ^Qj. assumed inclination angle of 



Emerge ^ 511_^4j 
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merge 



' 267!i,f Myr 



Taken at face value, the detection of this remarkable bi- 
nary implies that the number of such systems in the Galaxy 
is very large. The WD has an apparent g-band magnitude 
« 2 magnitudes brighter than the limiting magnitude of the 
survey, 18.9 (Badenes et al. 2009; Muflally et al. 2009). 
Thus, the volume probed by the survey for systems analogous 
to SDSS 1257+5428 is Vswarms - /sDSs(47r/3)(125pc)3 ~ 
2 X 10^ pc-', where /sdss ~ 1/4 is the fraction of the sky 
probed by SDSS to g « 18.9 in SWARMS, and where we have 
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taken /sDSs(47r/3) ^ 1. The total volume of the Galaxy is 
Vmw ~ Inr^h, where r ^ 8.5r8.5kpc kpc is the galacto-centric 
radius and h 4/i4kpc kpc is the scale height accessible to ob- 
jects with total space motion of ~ 100 km s"'. Thus, 



Vswarms 
Vmw 



1 X 10 *'"8.gkpc''4kpc- 



(1) 



The primary point of this Letter is that this ratio is exceed- 
ingly small. Since equation (1) is equivalent to the probability 
of detecting such a system if it was the only such system in the 
Galaxy, we conclude that it is not the only such system. More 
probable is that each of the VswARMS-sized volumes that con- 
stitute the Galaxy is populated by > 1 such binaries. The 
total number of WD-NS/BH binaries in the Galaxy is then*" 



and the current merger rate is 



I ^ iU-yV 5kpeft4kpc (2) 

V Vswarms / 
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flife 



.5xlO-%yr-'ri5kpc/^4kpcf2Gyr, (3) 



where ficyr = fiife/2Gyr and = fcooi+fmerge (e.g., Phinney 
1991). 

We caution that the estimate of the number of such systems 
in the Galaxy in equation (2) and the merger rate in equation 
(3) is highly uncertain. Ignoring the uncertainties in the ratio 
of volumes and in fyfe, the Poisson error on A^, given the single 
detection, yields a lower Hmit of N > 0.05 (A^mw > 5 x 10"*; 
Fmw > 2.5 x lO'^ yr"') and > 0.01 (A^mw > 1 x IO^Fmw > 
5 X 10"^ yr"') at 95% and 99% confidence, respectively (e.g., 
Gehrels 1986). A single additional such binary detected in the 

* We note that the density of stars in the Solar neighborhood can be esti- 
mated similarly, and with at least order-of-magnitude fidelity, from the first 
measured parallax of 61 Cygni by Bessel (1838). 
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local ^10*" pc^ would provide strong evidence for '--^ 1 . Un- 
fortunately, the lack of such a future detection in SWARMS 
does not strongly constrain Fmw-^ In addition to the uncer- 
tainty in N, there is uncertainty at the factor of ~ 2-4 level 
in Vmw; however, our fiducial values for r and h are reason- 
able, given the age and space velocity of the binary consid- 
ered. Note also that SWARMS selects for edge-on binaries 
due to the low resolution of the SDSS spectra used to tar- 
get interesting WDs via radial velocity variations (Badenes et 
al. 2009; MuUally et al. 2009). This effect may increase the 
inferred rate by a factor of order ~ 2 — 4. The ambiguity in 
the cooling age of SDSS 1257H-5428 adds an additional factor 
of ~ 2 - 4 uncertainty. Finally, our estimate tacitly assumes 
that the binary is equally detectable by SWARMS over its few 
Gyr lifetime from birth to death; the fact that the WD was 
brighter in the past implies a lower overall rate than given in 
equation (3), but the magnitude of this correction depends on 
the luminosity function of WDs in such binaries. 

2. DISCUSSION: THE RATE 

The nominal estimate of the merger rate for WD-NS/BH 
binaries in equation (3) is very high. Taking the core-collapse 
supernova (SN) rate to be ~ 1 - 2 x 10~^yr~^ in the Galaxy, 
we see that the WD-NS/BH merger rate is ~ 20 - 40 times 
smaller. Comparison with the observed rate of core-collapse 
SNe in the local 100 Mpc volume (see, e.g., the compilation 
of data by Horiuchi et al. 2009), ~ 1 - 2 x 10^ Gpc'^ yr"', 
implies that the average rate of WD-NS/BH mergers in the 
local universe is ^ 0.5 -1x10^ Gpc"^ yr"' . This is a factor of 
just ~ 3 - 6 lower than the local Type la rate, ~ 3 x 10'* Gpc"^ 
yr-i (Dilday et al. 2008). See Table 1. 

Although the rates of short- and long-duration gamma-ray 
bursts (GRBs) are uncertain because of the overall beam- 
ing correction for the two populations, Fmw in equation (3) 
is at least ~ 10^ times larger than the rates inferred for ei- 
ther population; the long-duration GRB rate is observed to be 
^ 1 Gpc"^ yr~', and the beaming-corrected rate is ^ 30Gpc"^ 
yr"' (e.g., Guetta et al. 2005). The rate of short-duration 
GRBs is also small (~ lOGpc'^ yr"'; Nakar 2007). This dis- 
parity between Fmw and the GRB rate is particularly inter- 
esting because these mergers have been suggested as a mech- 
anism for some long- and short-duration GRBs (§3). Inter- 
estingly, our rate for Fmw is only ^ 5 times higher than the 
current estimate for the rate of NS-NS mergers in the Galaxy 
(Kalogera et al. 2004ab), often thought to be associated with 
short-duration GRBs (see Table 1). 

Our estimate of Fmw is also significantly larger than the 
rate derived from statistical methods based on the 3 previ- 
ously known WD-NS(pulsar) binaries that will merge within 
a Hubble time, ^ 0.2- 10 x 10~^yr~', obtained without mak- 
ing an overall correction for pulsar beaming (Kim et al. 2004). 
This result is statistically dominated by the single WD-NS bi- 
nary system PSR Jl 141-6545, whose merger timescale is sim- 
ilar to the system described here. The rate derived in Kim et 
al. (2004) is about an order of magnitude smaller than that 
obtained from population synthesis studies (e.g., Kalogera et 

' Based on the constraints on the mass function alone, there is a signif- 
icant chance (~ 50%) that one of the systems recently reported by Kilic et 
al. (2009) is also a WD-NS binary. Since the survey volume is much larger 
for the Kilic et al. survey (they target bright Helium WDs), if the system is 
confirmed to be a WD-NS binary, Fmw is still likely to be dominated by 
SDSS 1257+5428. 



TABLE 1 
Events & Rates 



Event Type 


Rate 


Reference 




(10'*Gpc-3 yr-i) 




WD-NS/BH Merger 


0.5-1 


this work 


Type II SN 


10-20 


Horiuchi et al. (2009) 


Type la SN 


3 


Dilday et al. (2008) 


NS-NS Merger 


CP 


Kalogera et al. (2004ab) 


Short GRB 


0.001''''= 


Nakar (2007) 


Long GRB 


0.0001* 


Guetta et al. (2005) 



^Calculated using the rate fa 10"^ yr ' derived in Kalogera et al. (2004ab) 
and the conversion factor of 10"^ Mpc"' based on the B-band luminosity of 
the local universe from Kalogera et al. (2001), their Section 4. 

''Uncorrected for beaming. 

'^Short-duration GRBs with extended emission (e.g., 060614) amount to ~ 
25% of the total short-duration GRB population (Norris et al. 2009). 

al. 2005), and it is about ~ 10^ times smaller than Fmw in 
equation (3). The estimate by Fryer et al. (1999a) is similiirly 
small, but with large uncertainties. 

Our estimate of Fmw is ~ 1 - 10 times higher than the rate 
advocated by Davies et al. (2002) in their binary evolution 
scenario developed to explain the WD-pulsar systems Jl 141- 
6545 and B2303H-46.** In their model, the initial primary 
transfers enough of its mass to the secondary that the latter 
becomes the more massive star. The core of the original pri- 
mary makes the WD, which enters a common envelope with 
the more massive secondary when it evolves off the main se- 
quence. In the case of their PSR J1141-6545-like binaries, 
the secondary helium star evolves to contact with the primary 
and the majority of its envelope is ejected from the system. 
The secondary then undergoes core collapse, producing a NS. 
However, since SDSS 1257H-5428 is not highly eccentric it 
is unclear to what extent it may be viewed as analogous to 
Jl 141-6545 (see [10], §3). 

Our 95% confidence lower limit on the rate — Fmw ^ 
2.5 X 10"^yr~' — is still larger than either of the observed 
or beaming-corrected short- or long-duration GRB rates, but 
overlaps with the probable distribution of NS-NS merger rates 
for the Galaxy (Kalogera et al. 2004ab), the upper end of the 
pulsar beaming-corrected rates derived from PSR Jl 141-6545 
(Kim et al. 2004), the rates inferred from population synthesis 
studies (Davies et al. 2002; Kalogera et al. 2005), the rate de- 
rived by Nelemans et al. (2001), and the estimate of Edwards 
& Bailes (2001) of ~ 10"^ yr"' based on PSR Jl 141-6545. 

3. IMPLICATIONS 

The estimate of the merger rate in equation (3) has a number 
of implications: 

1. GRBs & Transients: The merger of WD-NS/BH bina- 
ries has been discussed as a mechanism for producing long- 
duration GRBs (e.g.. Fryer et al. 1999b) and the new class 

Note that their production rate of tight binaries is highest (closest to eq. 3) 
for a flat (i.e., not Salpeter) distribution of secondary masses and for large 
values of the common envelope evolution parameter. 
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of short-duration GRBs with extended emission (e.g., GRB 
060614; King et al. 2007; see Gehrels et al. 2006 & Gal- Yam 
et al. 2006). A recent study by Norris et al. (2009) finds that 
only ~ 25% of short duration GRBs have extended emission. 
Thus, if King et al. (2007) are correct and WD-NS mergers are 
the central engine for 060614-like bursts, then either equation 
(3) overestimates the rate by a factor of ^ 10- - 10-', only a 
very small percentage of WD-NS/BH mergers produce GRBs, 
or these events are very highly beamed — with jet opening 
angles many times smaller than most estimates (e.g., Nakar 
2007; Guetta et al. 2005). Similar statements can be made for 
short-duration GRBs without extended emission, or for long- 
duration GRBs unaccompanied by SNe (see Table 1). The 
natural timescale for such transients is long (> 100 s) since 
the WD will be tidaUy disrupted at a radial distance from the 
NS comparable to its own radius. Equation (3) implies that 
transients generated by WD-NS/BH mergers should be rela- 
tively common in Milky Way-hke galaxies (~ 10-30% of the 
la rate), contrary to what is observed for long-duration GRBs 
(Stanek et al. 2006). Whatever transients these mergers pro- 
duce (see [4] below), they should be associated with relatively 
old Gyr) stellar populations, may occur many kpc from 
their host galaxies, and should not correlate with recent star 
formation. 

2. Potential Neutrino Signature: Whether the unseen compan- 
ion in SDSS 1257H-5428 is a NS or BH, the merger may gen- 
erate an interesting neutrino signature, with a maximum pos- 
sible total energy radiated of ~ 3 x lO^"' ergs (similar to a NS- 
producing core-collapse SN). However, because the accretion 
rate is likely to be modest, < 10"^ M0 s~', the merger disk is 
unhkely to be radiatively efficient in neutrinos (e.g., Popham 
et al. 1998; Chen & Beloborodov 2007). Since BH formation 
is a likely outcome of the merger (see [6] below), a distinct 
signature of this NS-to-BH transformation may be seen in the 
neutrino (or potentially high-energy photon) Ughtcurve. 

3. Nucleosynthesis: The outflow produced from the NS/BH 
during the accretion phase of the WD — either in a jet or 
a wind — may produce thermodynamic conditions favorable 
for interesting nucleosynthesis. For very high accretion rates 
(- 1 Mq s"'; Chen & Beloborodov 2007; Metzgeret al. 2009) 
the outflow may become neutron rich, producing conditions 
favorable for production of the r-process nuclei (e.g., Met- 
zger et al. 2007). Even in the absence of a neutron excess 
(as is more likely for the low accretion rates one estimates for 
such a merger/disk) if the dynamical timescale for ejection is 
short enough, the r-process might still occur (Meyer 2002). 
Since the rate estimated in equation (3) is relatively high with 
respect to NS-NS mergers (a commonly discussed r-process 
production site; Freiburghaus et al. 1999), WD-NS/BH merg- 
ers may make an interesting contribution to the heavy-nucleus 
budget of the Galaxy. 

4. Optical Signature: Similar to the case of WD-WD merg- 
ers, or accretion-induced collapse of a WD to a NS as a re- 
sult of accretion from a companion, as considered in Metzger 
et al. (2009), there may be a low-luminosity optical transient 
generated by < 10"^ M0 of ^^Ni associated with the merger of 
WD-NS binaries. This material may be ejected as an initially 
proton-rich outflow from the disk around the NS or BH. The 
total ejecta mass would be small on the scale of normal SNe. 
The rate in equation (3) is small enough that such transients 
may be so far unknown, but we speculate that they might 
be Unked to classes of unusual SNe like SN 1991bg (Filip- 



penko et al. 1992), 2002bj (Poznanski et al. 2009), 2008ha 
(Foley et al. 2009; Valenti et al. 2009), or 2005cz (Kawa- 
bata et al. 2009) and 2005E (Perets et al. 2009). Such events 
should be relatively common in the local universe (^ 10-30% 
of the la rate), trace relatively old stellar populations, have 
no optically- visible progenitors, be seen many kpc from their 
host galaxies (e.g., SN 2005E), and many should be seen by 
current transient survey efforts such as LOSS (Li et al. 2000), 
Palomar Transient Factory (PTE) and the Panoramic Survey 
Telescope and Rapid Response System (PanSTARRS); the 
future Large Synoptic Survey Telescope (LSST) should see 
many hundreds per year. The composition of the ejecta may 
vary, depending on the type of WD disrupted (e.g.. He, C/O, 
O/Ne/Mg). 

5. Gravitational Waves: Assuming Mwd = 0.9 Mq, Mns/bh = 
1.6 Mq, and cosj = for SDSS 1257H-5428 results in a 
gravitational wave (GW) strain ampUtude of w 2.3 x 10"^', 
which exceeds the nominal sensitivity of LISA of ^ 10"^' 
at /gw ~ 0.1 mHz (Roelofs et al. 2007). The estimate of 
equation (3) suggests that WD-NS/BH binaries will also con- 
tribute significantly to the GW background in the LISA band 
(Badenes et al. 2009; Kim et al. 2004), which may affect the 
detectability of the SDSS 1257H-5428 signal. Importantly, in- 
dividual mergers in the local universe may also be important 
GW sources for LIGO if the merger results in BH formation 
or if the massive NS merger remnant is initially rapidly rotat- 
ing (Garcia-Berro et al. 2007; Paschalidis et al. 2009). 

6. The Lowest Mass Black Holes: Note that the total mass of 
SDSS 1257H-5428 probably exceeds the maximum mass for 
a NS (e.g., Lattimer & Prakash 2007). Thus, if flie WD was 
entirely accreted, this event would be accompanied by BH 
formation with ~ 2.5 M0 (see Brown et al. 2001). The rate of 
WD-NS/BH mergers estimated in equation (3) is large enough 
that the Galaxy should be littered with ~ 2-3 M© BHs; there 
should be ^ 10*- 10' low-mass free-floating BHs throughout 
the Galaxy. 

7. The Origin of Ultra-High Energy Cosmic Rays: If WD- 
NS/BH mergers produce a BH and accretion disk (see 
Popham et al. 1998; Fryer et al. 1999b), the magnetic lu- 
minosity via the Blandford-Znajek mechanism will hkely be 
large enough for production of UHECRs (e.g., Waxman & 
Loeb 2009; see also Waxman 1995). Because the energy 
reservoir in WD-NS mergers is comparable to the energy bud- 
get of GRBs, but the overall rate of WD-NS/BH mergers 
within the local GZK volume (~ 100 Mpc) is > 100 times 
larger (eq. 3; §2), WD-NS/BH mergers may dominate the 
UHECR budget. The merger may also produce a rapidly rotat- 
ing (ms spin period) NS with potentially short-lived magnetar- 
strength (^ lO'^ G) magnetic fields. In this case, the mecha- 
nism advocated by Arons (2003) for UHECR production may 
obtain. The rate needed by Arons to account for the observed 
UHECR budget is close to the rate derived for WD-NS/BH 
mergers estimated here (Table 1). Additionally, note that 
the formation of ms magnetar-like conditions via WD-NS/BH 
mergers alleviates the problem of getting UHECRs out of the 
overlying dense massive star progenitor discussed by Arons. 

8. Binaries in Type Ib/c SNe: If the scenario outlined by 
Davies et al. (2002), which predicts a formation rate of SDSS 
1257-H5428-like binaries within a factor of ~ 1 - 10 of the es- 
timate in equation (3), is correct, then the star that produces 
the NS (originally the secondary) explodes after becoming a 
He star and transferring a significant fraction of its envelope to 
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the WD, and potentially expelling it from the system. The fact 
that the overall rate of Type-Ib/c SNe is ^ 10-20% of core- 
collapse SNe (Prieto et al. 2008) implies that if the Davies et 
al. mechanism is correct, and if Fmw is correct, then many 
Type Ib/c SNe (~ 20-50%) explode with a very close WD 
companion. 

9. Tight Companion Interaction: Regardless of the binary for- 
mation channel, the SN explosion that produces the NS may 
interact with the close secondary, be it a WD, main sequence 
star, or otherwise. In this case, there may be a signal of in- 
teraction in the very early-time hghtcurve of many stripped 
envelope Type-Ib/c SNe as a result of the break-out flash and 
Shockwave interacting with the nearby companion (see, e.g.. 
Marietta & Burrows 2000; Kasen 2009). 

10. Star & Binary Formation: The Galactic birth rate of mas- 
sive stars is approximately equal to the core-collapse SN rate. 
Thus, equation (3) implies that ^ 2-4% of all massive stars 
are born with a close binary companion capable of produc- 
ing a SDSS 1257H-5428-like binary. If true, this has important 
impUcations for the physics of massive star formation and the 
demographics of young massive star binaries (e.g., Krumholz 
et al. 2009), as well as for pulsar binaries. The low eccen- 
tricity of the orbit of SDSS 1257H-5428 suggests interaction 
with the companion subsequent to the NS birth. Addition- 
ally, if WD production precedes the NS, equation (3) implies 
a significant number of young pulsars with very tight WD 
companions. However, only one such system has been de- 
tected among the > 1000 non-millisecond pulsars (MSPs) — 
the highly-eccentric PSR JI 141-6545 (see Edwards & Bailes 
2001; Davies et al. 2002; Kim et al. 2004) — disfavoring this 
as an analog. If the NS preceded the WD and was recycled 
into a MSP, then MSPs should be embedded in the centers 
of a fraction of young planetary nebulae, potentially visible 
either in radio and/or X-rays and gamma-rays. 



Indeed, the estimate of equation (3) is large enough that the 
simplest explanation for the discovery of SDSS 1257H-5428 is 
either that SWARMS was very lucky, or that instead of a WD- 
NS/BH binary, this system is a tight WD- WD binary (Badenes 
et al. 2009). Although the high nominal space velocity of the 
system argues against this possibility, if SDSS 1257H-5428 is 
in fact a WD-WD binary, then the similarity between the de- 
rived rate and the observed Type la SN rate is not a coinci- 
dence (see Table 1), particularly in light of the results of Mul- 
lally et al. (2009) (see also Kihc et al. 2009). The relatively 
high value for the rate in equation (3) is yet more perplexing 
if the companion to the observed WD is in fact a BH, and not 
a NS (see Brown et al. 2001). Nevertheless, the relative lack 
of pulsars with tight binary companions may argue for a BH 
companion, and it should be kept in mind that the SWARMS 
survey is perhaps the first where such tight WD-BH binaries 
could have been detected. Clearly, a more complete census 
and analysis — as will be provided by the SWARMS survey 
(Mullally et al. 2009) — and follow-up observations of SDSS 
1257-1-5428, including a search for (potentially pulsed) radio 
and X-ray emission as well as a parallax measurement, are 
needed. These will be vital, since, if SDSS 1257H-5428 is 
radio-quiet, this would imply a class of binaries distinct from 
any presently known. 
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